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SUMMARY 

An analysis of the steady-state aerodynamic heating problem 
at high-supersonic speeds is made for two-dimensional flows with 
laminar boundary layers. The aerodynamic heating is shown 
to be reduced substantially by injecting a small amount of coolant 
through a porous surface into the boundary layer. The coolant 
injection acts in two ways to decrease the aerodynamic heating: 
First, and most important, the velocity profile is altered such 
that the rate at which heat is conducted to the surface is reduced 
and, second, the coolant absorbs an amount of heat which is a 
function of the difference in temperature between the surface 
and the coolant. The first effect provides the advantage of cooling 
by injection over that of simply using a coolant to absorb heat 
from the surface. Calculations of the stability of the laminar 
boundary layer show that for a wide range of high-speed flight 
conditions the boundary layer would remain laminar at all 
Reynolds numbers according to the stability-theory consideration. 
The analysis includes calculations of the cooling requirements 
and equilibrium surface temperatures for flat plates and for 
flat porous surfaces with several rates of fluid injection at Mach 
numbers from 5 to 15 and altitudes from sea level to 200,000 
feet. Some calculations of the skin friction are also included. 

INTRODUCTION 

The aerodynamic heating problem assumes considerable 
importance at high-supersonic speeds. Sanger and Bredt 
(reference 1) have calculated the high-speed aerodynamic 
forces and equihbrium surface temperatures at extremely high 
altitudes where the molecular mean free path is large (free- 
molecule-flow region) compared with a characteristic body 
dimension. Although under these conditions the surface 
temperatures are low, the maximum lift-drag ratios become 
very small. Consequently, flight at- these altitudes maybe 
satisfactory for nonlifting missiles but would be unsatisfactory 
for steady level flight. On the other hand in the lower atmos- 
phere (say, sea level to 200,000 ft) the lift-drag ratios are 
considerably larger, but the aerodynamic heating is most- 
acute and some means of surface cooling must be employed 
at very high speeds. The scope of the present report is 
limited to a consideration of the steady-state aerodynamic 
heating problem at altitudes where the air may be considered 
a continuum, and particular emphasis is placed on a means of 
cooling the surface by continuous injection of a fluid in the 
boundary layer. 

Both the heat transfer and the drag coefficients are known 
to be lower for laminar than for turbulent flows. At super- 
sonic speeds it may be possible to maintain a laminar bound- 
ary layer and thus alleviate the aerodynamic heating 


somewhat. The theoretical investigation of Lees (refer- 
ence 2) on the stability of the laminar boundary layer in 
compressible flow indicates that the laminar layer is com- 
pletely stable at all Reynolds numbers at supersonic speeds 
for a sufficiently low ratio of surface temperature to stream 
temperature. Thus, for flow over smooth surfaces small 
disturbances are damped and a turbulent boundary layer 
does not develop even at high Reynolds numbers. The 
possibility of maintaining a laminar flow under these condi- 
tions with finite but small disturbances in the boundary 
layer — that is, maintaining a laminar boundary layer over 
surfaces that may be employed on a high-speed aircraft — has 
not been verified experimentally. Experimental verification 
of this question as well as the general conclusions of Lees* 
work is therefore desirable. Nevertheless, in view of the 
conclusion of laminar stability for infinitesimal disturbances, 
it is desirable to investigate theoretically methods of decreas- 
ing the aerodynamic heating for laminar flows. 

The aerodynamic heating is alleviated somewhat with 
increasing altitude, but even with laminar flow at altitudes 
of 200,000 feet the surface temperatures may be excessive at 
high-supersonic speeds. An investigation of means of cooling 
the surface is therefore desirable. One possible method of 
cooling is that of injection of a cool fluid through a porous 
surface into the boundary layer (reference 3). This coolant 
injection acts in two ways to decrease the aerodynamic 
heating: First, it alters the velocity profile such that the 
heat-transfer rate from the fluid to the surface is reduced 
and, second, the coolant absorbs an amount of heat which is 
a function of the difference in temperature between the 
surface and the coolant. The first effect provides the 
advantage of cooling by fluid injection over that of simply 
using a coolant to absorb heat from the surface. 

Injection of the coolant through a porous surface affects the 
boundary-layer stability in two ways. The direct effect of 
injecting a fluid in the boundary layer is to alter the velocity 
profile such that the flow is less stable. The indirect effect 
is to give a lower surface temperature which in turn tends 
to make the flow more stable. Actual calculations must 
determine whether the flow' would be stable for given 
conditions. This question will be discussed in more detail. 
The blowing rate should be kept as low as possible, consistent 
with adequate cooling, since the coolant will probably be 
carried by the aircraft. Consideration of storage would 
probably dictate the use of a liquid, and a coolant which has 
a low temperature and a high heat capacity would be most 
efficient. 
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In the present report, calculations of the steady-state 
heating of laminar flows are presented for several altitudes 
from sea level to 20.0,000 feet; Mach numbers of .5, 10, and 
15; and several different surface lengths. Results are pre- 
sented for a flat porous surface with three rates of coolant 
injection as well as for a flat plate (no blowing). For the 
case of continuous injection (injection through a porous sur- 
face) the calculations were based on a theory for Prandtl 
number of 1.0, and in all cases the results apply to a flat 
surface at zero incidence. Included in the present report 
are calculations of the stability of the laminar layer and 
some calculations of the skin friction. 


ANALYSIS 

In order to determine the equilibrium surface temperature 
of a body in motion through a fluid, a balance between the 
heat brought to and that removed from the surface must be 
made. This transfer of energy occurs mainly through the 
processes' of convection and radiation. The analysis is 
simplified by considering only nighttime operation. Thus 
the effects of solar and atmospheric radiation are neglected 
and only radiation from the aircraft surface is considered. 
With these restrictions, then, the heat balance may be written 
as 

H e +H r +AH=0 (1) 
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SYMBOLS 

average heat transfer rate at surface, 

Bt,u/(sq ft) (sec) 

rate at which heat is radiated from surface, 
Btu/(sq ft) (sec) 

rate at which heat is absorbed by coolant, 
Btu/(sq ft) (sec) 

heat rate to vaporize coolant, Btu/(sq ft) (sec) 
heat-removal rate, Btu/(sq ft) (sec) 
temperature, °F absolute 
length of surface from leading edge, feet 
coordinate norrrfal to surface, feet 

conductivity, Btu/(sec) (ft) (°F abs) 
heat capacity at constant pressure, 

Btu/ (slug) (°F abs) 
coefficient of viscosity, slugs/ (ft) (sec) 
mass density, slugs per cubic foot 
latent heat of vaporization, Btu/slug 
Stefan-Boltzmann constant, 

(4.8X10 _13 Btu/(sec)(sq ft)(°F abs) 1 ) 
velocity parallel to surface, feet per second 
velocity normal to surface, feet per second 


ratio of heat capacities 
stream Mach number 

Reynolds number 
Prandtl number 

coolant parameter 2 Pi u i' x j 

nondimensional stream function 
skin friction drag, pounds 
skin-friction drag coefficient 

value at surface 
coolant value 
stream value 
isothermal value 
stagnation value 
recovery value 


where H e is the average heat-transfer rate between the fluid 
and the surface, H r is the rate at which energy is radiated 
from the surface, and A H is the rate at which heal is ab- 
sorbed .by the surface. The temperature at which A//— 0 
is termed the natural equilibrium temperature. The radiant 
energy for a grey body is readily calculated for a given sur- 
face temperature T, (°F abs) from the expression 

B‘ r =0.9(rT, 4 (2) 

■where a is the Stefan-Boltzmann constant and the factor 0.9 
is chosen as the emissivity for the present calculations. The 
rate of heat radiation is given as a function of surface tem- 
perature in figure 1. The. determination ot the conductive 
heat rate requires the solution of the bouudary-layer equa- 
tions. The numerical solutions for high-speed flows are 
lengthy and tedious, and in order to expedite the work ap- 
proximate solutions were used. Although this procedure 
affects the absolute accuracy somewhat, the comparison be- 
tween the flat plate and a flat surface with fluid injection is 
unaffected. The solutions for the flat plate and the surface 
with coolant injection differ in some respects and are, there- 
fore, discussed separately. Specifically, the steady-stale 
heat transfer at supersonic speeds for a flat plate and for a 
flat porous surface- with fluid injection, both at a constant 
temperature over the entire surface, is considered herein. 

Conductive heat transfer.— Throughout the analysis the 
flow is considered laminar, and the temperatures thus calcu- 
lated are later compared with the results of the stability 
theory of the laminar boundary layer to determine their 
validity. The average heat rate between the fluid and a sur- 
face of unit width and of length x may be written according 



Figure 1.— Rata of boat radiation for emissivity at 0.8. 
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to Fourier’s law of heat transfer as 



where k is the conductivity of the fluid, T is the temperature, 
y is the coordinate normal to the surface, and the subscript s 
denotes the value at the surface. The quantity bT/by is 
found by solving the dynamic and energy equations of the 
boundary lat'er. Even for a simple co nfig uration such as a 
flat plate the problem as stated becomes very difficult for 
high-speed gas flows. Several investigators have calculated 
velocity and temperature profiles from the boundary-layer 
equations by various numerical integration methods(reference 
4, for example) ; however, these calculations are insufficient 
for the present analysis. The usual method of determining 
the high-speed heat transfer (references 5 and 6) is to modify 
the low-speed solutions. This practice is followed herein. 

Conductive heat transfer for a flat plate. — The usual gener- 
alization of the low-speed heat-transfer equations results by 
employing values of density, viscosity, conductivity, and heat 
capacity based on surface conditions instead of stream values 
and by modifying the temperature potential (reference 5). 
For laminar flow over a flat plate then the average heat- 
transfer rate between “the fluid and surface is 

H.= 0.66 '\pry*(T t -T rt ) (4) 

where the subscript 0 denotes the value for a flat plate, 
p, is the density at the surface, jtt, is the viscosity at the sur- 
face, Pr is the Prandtl number (taken as 0.72) , pi is the stream 
velocity, and T re is the recovery (or adiabatic wall) tempera- 
ture, that is, the temperature an insulated nonradiating sur- 
face would assume in an adiabatic flow. This form of the 
heat-transfer equation has proved satisfactory at moderate- 
supersonic speeds; however, it has not been verified experi- 
mentally at the high speeds considered herein. 

The ratio of the temperature rise with no heat transfer to 
the stagnation-temperature rise is termed the recovery factor. 
A knowledge of this factor permits the determination of the 
recovery temperature. In general this factor is dependent 
upon the geometry of the body and the Prandtl number and 
is found from experiment or calculation. A value of 0.9 is 
used herein and the recovery temperature then is found from 
the relation 

T„=T 1 +o.9(r,-:r 1 ) ( 5 ) 

where T x is the stream temperature and T t is the stagnation 
temperature. In order to allow for the variation in the ratio 
of specific heats with temperature, the stagnation tempera- 
ture was calculated from a relation derived by quantum 
statistical methods (reference 6) 

M 2 +5 _ T t 1.13 0.2 . 

6060 1212 ' r6060”l , r4050"l . w 

< “ p LTr_r 1 

where M is the stream Mach number. Figures 2 and 3 present 
the variation of stagnation and recovery temperatures 
respectively, with. Mach number for several altitudes. 
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The heat balance for laminar flow over a flat plate at 
several altitudes and Mach numbers was determined from 
equation (1). The altitude tables of reference 7 were used 
to determine the stream temperature and density, and the 
viscosity and conductivity coefficients were calculated from 
Sutherland’s equation. The heat absorbed by the surface 
(AH) for a 10-foot surface length is given as a function of the 
surface temperature in figure 4 for several altitudes. The 
natural equilibrium temperature is found from the condition 
that AH= 0. This temperature, for the flat plate (f 0 =0), 
is given in figure 5 for several altitudes, Mach numbers, and 
surface lengths. In general the equilibrium temperature is 
found graphically or by trial and error. 

Heat balance with continuous fluid injection. — The high 
natural equilibrium temperatures of a flat plate at high- 
supersonic speeds emphasize the need for some adequate 
means to cool the surface. Injection of a fluid through a 
porous surface into the boundary layer changes the velocity 


profile such that the velocity gradient 


du 

by 


at the surface is 


decreased. Since the temperature gradient, and thus the 
heat-transfer rate, is approximately proportional to the veloc- 
ity gradient, the rate at which heat is conducted to the surface 
is reduced. A very small rate of fluid injection decreases the 
aerodynamic heating significantly. Since the quantity of 
fluid injected is small, the usual boundary-layer equations 
apply. 

If the initial temperature of the fluid injected into the 
boundary layer is lower than the surface temperature, the 
coolant will absorb an amount of heat proportional to the 
temperature difference; furthermore", if the coolant is a 
liquid, it will absorb an additional amount of heat propor- 
tional to the heat of vaporization. The change in the 
velocity profile also alters the stability. of the boundary layer. 
In reference 3 it is shown that the fluid injection is destabi- 
lizing, but on the other hand the effect of the hep,t transfer 
is stabilizing. Thus, the net effect of injection and cooling 
must be accounted for in the stability consideration. For a 
wide range of conditions the flow is stable in the laminar 
form. See section entitled “Stability considerations of the 
laminar boundary layer.” 

For the present analysis it is assumed that the coolant is 
liquid air and that part of the heat in the boundaiy layer is 
utilized to vaporize the coolant. The heat balance then 
may be written as 


H c -\-H a -)-H c -\-H r =0 (7) 


The heat radiated from the porous surface is given by equa- 
tion (2) and the average value of H, is 
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where L is the latent heat of vaporization (taken as 2,834 
Btu/slug), E is the Reynolds number based on stream condi- 
tion, p, is the density of the coolant, v, is the velocity normal 
to the surface, and f 0 is a nondimensional flow parameter 
defined as 

(9) 

PittiV Ml 


The parameter f 0 arises in the treatment of the dynamic and 
energy equations (reference 3). The condition of constant 
temperature over the surface is satisfied for Constant. 

In a" similar manner the average rate at- which heat is 
absorbed by the coolant is found as 


J Q P>v,c Pi {.T,— T,) dx 
or 

H a =-r 0 ^|i c?i (T.-r,) do) 


where c fl is the specific heat and T t is the temperature of 
the vaporized coolant. The heat capacity a> t was taken as 
7.712 Btu/(slug) (° F abs) and the temperature T t , as 147° F 
absolute. 

As in the case of the flat plate, the calculation of the heat 
transfer between the boundary layer and the surface requires 
a solution of the bound ary-layer equations. However, for 
the case of continuous injection a solution with certain 
simplifying assumptions is readily determined from known 
isothermal solutions by a simple quadrature (reference 3). 


Then the quantity 



may be evaluated in terms of the 


slope of the velocity profile for the isothermal problem, and 
the average heat-transfer rate between the boundaiy layer 
and surface is readily found from equation (3). 

With the assumption that the Prandtl number is unity, 
the dynamic and energy equations have the same integral in 
the absence of a pressure gradient and the temperature is 
related to the velocity by 


T p, 

hi 

Tr p ~ 

Ti [ 
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( 11 ) 


where H a is the rate at which heat is absorbed by the coolant 
and H, is the heat rate required to vaporize the coolant. 


where w is defined as ujui and 7 is the ratio of specific heats. 



Heat-removal rate, A//, Btu /(sq ftjfsecj 
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0 fOO 200 300 400 500 600 700 800 900 tOOO O 100 200 300 400 500 600 700 800 900 1000 
Surface temperature, Z, °F Surface temperature, Ts, °F 



Surface temperature, T s , °F 


(a) Sea level. (b) Altitude, 50,000 feet. 

(c> Altitude, 100,000 feet. (d; Altitude, 150,000 feet, 

(e.i Altitude, 200,000 feet. 


Fiscal 4.— Average heat-removalrate for Sat plate of 10-foot length. 



Equilibrium surface temperature, T s , °F Equilibrium surface temperature. 
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Figure 5. — Equilibrium surface temperature for various amounts of blowing. 
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Figure 5 .— ‘Continued. 
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(e) Altitude, 200,000 feet. 
Figuse A— Concluded. 


Reference 3 stows that, for flow over a porous surface of 
uniform temperature with continuous injection, the equations 
of the compressible flow are of the same form as the iso- 
thermal case provided the Prandtl number is unity and the 
viscosity varies linearly with temperature. The coordinate 
n then is related to the isothermal value by 

+z=±a**)] (r-ro ft 

( 12 ) 


where the subscript iso refers to the isothermal value and 
f is a nondimensional stream function. From equations (11) 

and (12) the quantity (the value at w=0) can be ex- 

pressed in terms of the isothermal values of the slope of the 
velocity profile, and the average convective heat rate is then 
found as 



or 


H e =(H e ) to=0 



(13b) 


A quantity equal to four times 



is tabulated forseveral 


values of f Q in reference 8. The quantity (H’ e )f 0 =o is the 
convective heat rate for no blowing — that is, for a flat plate. 
This relation, however, does not give exact agreement with 
the results previously given for the flat plate. In order to 
provide a better basis for comparison of the equilibrium tem- 
peratures with and without blowing, the quantity (flc)r 0 =o 
may be replaced by H Ca ‘, therefore, the convective heat rate 
with blowing may be written as 



The heat balance for surfaces with coolant injection is 
determined from equation (7) by employing the foregoing 
relations. The equilibrium temperatures for several blowing 
rates, altitudes, and Mach numbers are given as a function 
of surface length in figure 5. 

Throughout this report, the boundary layer is assumed to 
be laminar, and this condition may not be fulfilled for some 
temperatures calculated. The results should then be checked 
against the stability criterions for laminar boundary layers. 

Stability considerations of the laminar boundary layer. — 
The direct effect of injection of a fluid in the boundary layer 
is to destabilize the flow; whereas, the effect of removing heat 
is stabilizing. The stability problem then is to determine the 
minimum critical Reynolds number (that is the Reynolds 
number at which the flow first becomes unstable) for a given 
surface temperature, Mach number, and amount of coolant 
injection. The range in which the heat-transfer calculations 
are valid can therefore be determined by comparing them 
with the results of such stability calculations. In reference 2, 
a detailed consideration of the stability problem is presented 
for compressible viscous fluids, and in reference 3, an analysis 
is presented for the specific case of flows over a porous surface 
with continuous injection. 

Instead of determining the minimum critical Reynolds 
number, it is more convenient to determine the surface 
temperature required to make the flow stable at all Reynolds 
numbers for a given Mach number and fluid injection rate. 
Because of the extreme curvature of the velocity profiles with 
injection of gas at the wall, the approximate formulas of 
reference 2 for the minimum critical Reynolds number are 
not applicable and the iteration method of that report must 
be employed. An estimate of the value of minimum critical 
Reynolds number, which will serve as a stability criterion, 
is obtained by taking the phase velocity to have its maximum 
possible value. This condition occurs when $<( 2 ) =0.580. 
(See reference 2 for the definition of the stability parameters.) 
For that value of $ 4 ( 2 ), the two branches of the curve of 
wave number a against R coincide and hence the phase 
velocity for each branch approaches the same value 




as I?-» 
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The actual calculation consists of choosing a stream Mach 
number and fluid injection rate. The value of c is taken as 

1 — -jg and a trial value of T s jTi is assumed. Equations (24) 

and (25) ofreference 2 are used to obtain a zeroth approxima- 
tion to the values of $<( 2 ) and $,( 2 ). Figure 9 and equations 
(28) and (29) of reference 2 are used to determine successive 
approximations to $<( 2 ). If the final value of <£,( 2 ) ^0.580, 
the assumed value of TJTi was incorrect and a new value 
must be tried with the corresponding new velocity profile. 

The surface temperature for stability is given as a function 
of Mach number for a given amount ol injection in figure. 6, 
and those., temperatures may be compared directly with 
those found from the heat-transfer calculations. If the 
calculated surface temperatures are equal to or less than 
those found from the stability considerations for a given 
Mach number and coolant injection rate, the flow is stable 
for all Reynolds numbers. For finite values of the minimum 
critical Reynolds number somewhat higher surface temper- 
atures for stability would be allowed. 

Skin friction. — The skin friction is. readily determined 
once the velocity profile is known. The skin friction drag is 
given by 



or, in terms of drag coefficient, 

r — 7? 4 ^ ps ( dw\ 

° f \ Pl u?x 


(15) 


For a linear variation of viscosity with temperature, equa- 
tion (15) reduces to 




Figure 8.— Critical surface temperature for boundary layers that are stable for all Reynolds 

numbers. 


The skin-friction drag coefficient calculated from equation 
(16) is presented in figure 7. This quantity is independent 
of Mach number because of the linear viscosity-temperature 
relationship. The skin-friction coefficient decreases with 
Mach number (Prandtl number less than 1.0) for a viscosity- 

/T \ 0 - 76 

temperature relation given by — =( jr) • (Soo r eference4.) 

The actual variation with Mach number, however, depends 
to a large extent on the nature of the viscosity-temperature 
relation. In general the variation of the drag coefficient 
with Mach number is small compared with the change due 
to fluid injection. 

DISCUSSION 

The results of the steady-state heating analysis for laminar 
boundary layers are shown in figures 4 and 5. Although the 
analysis refers specifically to a flat surface at zero incidence, 
the calculations are indicative of the high-speed heating 
problem in general. The calculations show that the steady- 
state temperatures are significantly dependent on the Mach 
number) and altitude and to a lessor extent on the surface 
length, and that continuous injection of a cool fluid in the 
boundary layer reduces the equilibrium temperatures 
substantially. 

The average amount of heat that must be removed from 
the surface of a flat plate 10 feet long to maintain a given 
equilibrium surface temperature is shown in figure 4. It is 
readily apparent that the cooling requirements are increased 
with high speeds and low altitudes. Similar calculations for 
other surface lengths indicate that the aerodynamic heating 
is alleviated somewhat for longer surfaces. 

The calculated equilibrium surface temperatures are pre- 
sented in figure 5 as a function of the length of the surface for 
several altitudes, Mach numbers, ami coolant injection rates. 
It "is apparent that the temperatures of the surfaces with 
coolant injection are appreciably lower than those of lire 
fiat plate (f 0 =0). For these calculations the coolant was 
assumed to be liquid air. By choosing a coolant with a higher 
heat capacity, even lower temperatures would result. From 
consideration of a high heat capacity and a low vapor tem- 
perature, liquid hydrogen would be a logical choice for the 
coolant. The injection of foreign gas may alter the stability 
considerations somewhat, and for that reason no calculations 
for the case in which hydrogen is used as the coolant are 
included. Figures 4 and 5 indicate that the equilibrium 
temperatures increase rapidly with Macli number and are 
highest for low altitudes. The effect of surface length on the 
equilibrium temperatures is not very significant. At the 
lower values of the coolant parameter f 0 the short lengths 
are somewhat hotter, but at the larger values of f 0 there is 
little variation with length. The amount of coolant injec- 
tion required to provide a sizeable decrease in the convective 
heat rate is relatively small. For example, for the most 
extreme, condition considered here (10-ft surface length, sea 
level, M= 15, |' 0 =— 1), the average mass flow is 0.00122 
slug/ (sq ft) (sec). At an altitude of 100,000 feet this value 
reduces to 0.00012 slug/(sq ft) (see). At a Mark number 
of 5 the mass flow is 0.0007 slug/(sq ft) (sec) and 0.00007 
slug/(sq ft) (sec) at altitudes of sea level and 100,000 feel. 
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respectively. For a flow parameter of f 0 = — 0-5 the mass flow 
would be one-half of those given. The average coolant mass 
flow decreases for increased surface lengths. 

The primary effect of coolant injection is the reduction in 
the amount of heat transferred from the boundary layer to 
the surface. This effect may be seen from equation (14). 
The heat-transfer rate with injection is equal to that for the 
flat plate times the ratio of the slopes of the velocity profiles, 
and this ratio is 0.492, 0.282, and 0.107 for f 0 equal to —0.5, 
—0.75, and —1.0, respectively. Thus, the heat rates at the 
surface for f 0 =— 0.5 and —1.0 are approximately }{ and Ko, 
respectively, of that of the flat plate.' 


The fact that the heat-transfer rate with injection is 
appreciably lower than that of the flat plate suggests the 
possibility that part of the boundary layer near the trailing 
edge could be sucked off and used as the fluid to be injected 
in the boundary layer. This air, of course, is hot and would 
heat up the boundary layer, but it might be expected that 
the net effect would be a reduction in surface temperature 
over that of a flat plate. A few calculations were made on 
this basis, but for all conditions calculated the heat added 
was found to be large and the net effect was to produce 
equilibrium surface temperatures slightly higher than those 
for a flat plate. 
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The critical ratio of surface temperature to stream temper- 
ature for stability of the boundary layer for all Reynolds 
numbers is given as a function of the stream Mach number in 
figure 6. Accurate calculation of the stability at high Mach 
numbers is difficult and for that reason the results are shown 
with dashed curves above a Mach number of 5. For all the 
results in figure 5, however, the stability analysis indicates 
that the boundary layer is laminar for all Reynolds numbers. 
The effect of injection on the stability is clearly shown in 
figure 6. At high-supersonic speeds (say Af> 3) the surface 
radiation alone provides adequate cooling to maintain a 
laminar boundary layer over smooth surfaces. 

The variation of skin-friction drag coefficient with Reyn- 
olds number is shown in figure 7. The fact that the skin 
friction is independent of Mach number arises from the 
assumption of a linear viscosity-temperature relationship. 
For other viscosity-temperature relations the. skin friction 
shows some variation with Mach number, and the actual 
variation depends on the nature of the viscosity-temperature 
relation. The skin friction is. relatively small and decreases 
with increased blowing rates as long as the boundary layer is 
laminar. 

CONCLUDING REMARKS 

An analysis of the steady-state aerodynamic heating 
problem at high-supersonic speeds for two-dimensional flows 
with laminar boundary layers shows that the aerodynamic 
heating is reduced substantially by injecting a small amount 
of coolant through a porous surface into the boundary layer. 
Calculations of the stability of the laminar layer show that 
for moderate rates of injection the boundary layer is stable 
at all Reynolds numbers, and these injection rates provide 
adequate cooling for many high-speed flight conditions. The 
calculations show that the equilibrium surface temperatures 
are significantly dependent on the Mach number and altitude 
and to a lesser extent on the surface length. In view of the 


conclusion of laminar stability and the large decrease in aero- 
dynamic heating with a small amount of coolant injection, 
it is desirable to investigate the problem further both experi- 
mentally' and analytically since the method appears practical 
for some high-speed flight applications. Further investiga- 
tions should include a study of the effectiveness of injecting a 
fluid into the laminar sublayer to decrease turbulent aero- 
dynamic heating. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., September 29, 19/ f 9. 
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